// Accomminotti and Flandreau 2008 WP

// Load original dataset
clear
insheet using "AF2008 WP Rep Data.csv"

//REPLICATION
summarize
display c(k)
mdesc
tabmiss
codebook codei
codebook year
gen gdppci = gdpi/popi
gen gdppcj = gdpj/popj
gen gdppcavg = (gdppci + gdppcj)/2
summarize

// Table 4
// Model 1
reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij

// Model 2
reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout

// Model 3
reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout an1 an2 an3 an4 an5 an6
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout an1 an2 an3 an4 an5 an6

// Model 4
reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij ymfnij an1 an2 an3 an4 an5 an6
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij ymfnij an1 an2 an3 an4 an5 an6

// Model 5
xtset codeij year
xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij, re
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij

// Model 6
xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout, re
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout

// Model 7
xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout an1 an2 an3 an4 an5 an6, re
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout an1 an2 an3 an4 an5 an6

// Model 8
xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij ymfnij an1 an2 an3 an4 an5 an6, re
codebook year if e(sample)
codebook codei if e(sample)
mdesc ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij ymfnij an1 an2 an3 an4 an5 an6

// IMPUTATION
clear
use "AF2008 WP Imp Data.dta"
summarize

mi import flong, m(imp) id(year codeij) imp(codei-codej impyeari-ln_distij)

mi xtset codeij year

// Table 4
// Model 1
mi estimate, esampvaryok post: reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij

// Model 2
mi estimate, esampvaryok post: reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout

// Model 3
mi estimate, esampvaryok post: reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout an1 an2 an3 an4 an5 an6

// Model 4
mi estimate, esampvaryok post: reg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij ymfnij an1 an2 an3 an4 an5 an6

// Model 5
mi estimate, esampvaryok post: xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij, re

// Model 6
mi estimate, esampvaryok post: xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout, re

// Model 7
mi estimate, esampvaryok post: xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij mfnij mfnout an1 an2 an3 an4 an5 an6, re

// Model 8
mi estimate, esampvaryok post: xtreg ln_impij ln_gdpi_gdpj ln_gdpij_popij ln_distij borderij langij commonij itaij ausij ymfnij an1 an2 an3 an4 an5 an6, re
